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Random walks on weighted networks 
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Random walks constitute a fundamental mechanism for a large set of dynamics taking place on 
networks. In this article, we study random walks on weighted networks with an arbitrary degree 
distribution, where the weight of an edge between two nodes has a tunable parameter. By using 
the spectral graph theory, we derive analytical expressions for the stationary distribution, mean 
first-passage time (MFPT), average trapping time (ATT), and lower bound of the ATT, which is 
defined as the average MFPT to a given node over every starting point chosen from the stationary 
distribution. All these results depend on the weight parameter, indicating a significant role of 
network weights on random walks. For the case of uncorrelated networks, we provide explicit 
formulas for the stationary distribution as well as ATT. Particularly, for uncorrelated scale-free 
networks, when the target is placed on a node with the highest degree, we show that ATT can 
display various scalings of network size, depending also on the same parameter. Our findings could 
pave a way to delicately controlling random-walk dynamics on complex networks. 

PACS numbers: 05.40.Fb, 89. 75. He, 05.60.Cd, 05.40.-a 



I. INTRODUCTION 



Random walks on complex networks have attracted 
increasing interest in the past few years P, Q because 
of their wide applications in a large variety of scientific 
fields j3j. Thus far, random walks on binary networks 
have been extensively studied, leading to important un- 
derstanding of this paradigmatic dynamical process and 
uncovering the effects of the underlying network struc- 
ture on the behavior of random- walk dynamics [4] . For 
example, general solutions for mean first-passage time 
(MFPT) from one node to another in a general network 
have been found by using different methods [E'-?]. The 
MFPT from node i to node j, denoted by Fij, is defined 
as the expected time for a walker starting from node i 
to first reach node j. This quantity is required in vari- 
ous considerations [Sj. Recently, average trapping time 
(ATT), (Fj), for the trapping problem a particular 
random walk defined as the average of Fij over all start- 
ing points i to a given trap node j, has received consider- 
able attention [10|-|28| , for this quantity can be applied to 
characterize transport efficiency. There are different ap- 
proaches for evaluating (Fj), showing nontrivial (depen- 
dence) scalings of MFPT {Fj) with network sizes, which 
strongly depend on network structural properties such 
as power-law degree distribution [H, [13: fractal- 
ity [lil, and modularity [13, [H- 

It is well known that in addition to the above- 
mentioned important structural properties, most real- 
life networks exhibit a large heterogeneity in the distri- 
bution of edge weights, characterizing the intensive in- 



teractions between node pairs |29|, |30[. The inhomoge- 
neous weights have a substantial effect on various dy- 
namical processes running on a network, including zero 
range processes [Slj, synchronization_ |32|. Issl . cascad- 
ing failures 



34| , mass-aggregation 
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^ ^, ^ , ^, diffusive pro- 

cesses [36| , bimolecular chemical reactions [37*1 , diffusion- 
annihilation processes [S^ , voter behavior ^391] , and traf- 
fic fluctuation ji^]. Thus far, a comprehensive analysis 
about random walks of weighted networks is still lacking, 
although it is suggested that weight heterogeneity could 
play an essential role in random-walk dynamics. 



In this paper, we present a rather comprehensive study 
of random walks on a generic weighted network, where 
the weight Wij of the edge linking nodes i and j has strong 
correlations with the degrees, di and dj, of nodes i and j, 
and is assumed to be Wij = (didj)^ with 9 being a control- 
lable parameter. Such a form of edge weights can be ob- 
served in various real weighted networks [29', . Based 
on spectral graph theory |41|, we will derive analytically 
the formulas of the stationary distribution, MFPT Fij 
from node i to node j, ATT {Fj) to a given target node 
j for a walker starting from a source node selected from 
the stationary distribution, as well as a lower bound for 
{Fj). The obtained results show that all these quantities 
are strongly affected by network weights characterized 
by parameter 9, which indicates the significant impact of 
weights on random walks. For the special case of uncorre- 
lated networks [i^l, we will provide succinct expressions 
for some interesting quantities. Particularly, for the trap- 
ping problem on an uncorrelated scale- free network, when 
the perfect trap is placed at a hub node with the largest 
degree, by tuning the weight parameter 9, ATT can ex- 
hibit different scalings with network sizes. Our work thus 
might have practical implications for controlling random 
walks on real- world networks. 
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II. FORMULATION OF THE PROBLEM 

Most existing works on random walks are confined to 
binary networks For a connected binary network G 
with N nodes and M edges, where the nodes are labeled 
by 1, 2, 3, . . . , A^, its connectivity is represented by the 
adjacency matrix A, in which the entry aij = 1 (or 0) if 
nodes i and j are (not) adjacent. The degree of node i 
is di = X^iLi : ^-iid the diagonal degree matrix D of 
network G is defined as follows: the «th diagonal entry 
is di, while all the nondiagonal entries are equal to zero. 
It is easy to verify that the total degree of all N nodes 
IS. K = 2M — X^jLi d'ij the average node degree is 
(d) = 2M/N. 

For discrete-time random walks defined on a binary 
network, at each time step the particle (walker) starting 
from its current position moves to each of its neighbor- 
ing nodes with the same probability. Such a stochas- 
tic process is characterized by the transition matrix 
P — D~^A |6|], in which the entry pij = aij/di presents 
the jumping probability from node i to node j in one 
step. This process can be described by an ergodic 
Markov chain [i^ |3] ; whose stationary distribution is 
TT = (tti, 7r2, . . . , ttac)^, where tt^ = di/K. It is clear that 
TT^P = TT^ and X^ili = 1- 

The above-defined random walks can be easily ex- 
tended to be on weighted networks [i^ . Mathematically, 
topological and weighted properties of a weighted net- 
work are described by a generalized adjacency matrix W, 
whose entry Wij specifies the weight of the edge connect- 
ing nodes i and j. In this paper, we focus on an undi- 
rected network having symmetric nonnegative weights 
Wij = Wji > 0. Moreover, we assume the weight of the 
edge between nodes i and j to be Wij = [didjY , where 9 
is a controllable parameter. Our assumption is based on 
empirical work on realistic networks, including scientific 
collaboration networks (29|] , metabolic networks [s^l , and 
airport networks [1^ l3Cll|. Moreover, the strength Si of 
node i is defined by Si = X^jen "^ij pSj. where the sum 
runs over the set fi^ of the neighboring nodes of i. 

Notice that random walks occurring on weighted net- 
works are biased. In the process of a random walk, the 
transition probability pij from node i to j is given by 



j / Si, which constitutes an entry of the transition 



matrix P = S'^^W for the biased random walk, where 
iS" is the diagonal strength matrix with its ith diagonal 
entry equal to the strength Si of node i. In the following, 
let s denote the sum of strengths for all the N nodes, 

i.e., s = X^ili^i = Sill SiLi ^ij- It is clear that the 
transition matrix P is a stochastic matrix. 

Next, we will show that many quantities related to the 
biased random walks are encoded in the transition matrix 
P and thus can be derived from P. Generally, except 
regular networks, P is asymmetric. So, we introduce the 
matrix 



which is real and similar to P and thus has the same set 
of eigenvalues as P. Furthermore, if is an eigenvector 
of matrix F associated with eigenvalue A, then S'^^ip is 
an eigenvector of P corresponding to eigenvalue A. 

It is easy to verify that for any given node i, J2]=i Pij — 
1 always holds. Thus, 1 is an eigenvalue of transition 
matrix P, as well as matrix F. Moreover, it is the great- 
est eigenvalue with single degeneracy. For eigenvalue 1 
with a corresponding normalized eigenvector iv, one has 
Tw = w, which together with Eq. ([1]) gives P{S~^w) = 
{S^^w). Because PI = 1, one has S~^w = 1. Thus, w 
can be determined explicitly as 

w = [wi,W2,-..,wn) = ( V ~' V ~' ■ • ■ ' y — ) • 

(2) 

Equation ([I} shows that matrix F is real and symmet- 
ric; thus all its eigenvalues are real. Let Ai, A2, A3, • • • , 
Xn be the N eigenvalues of matrix F for a network of size 
N, rearranged asl = Ai>A2>A3>...>Ajv>— 1, 
and let ipi, -02, "03 1 • • -j i'N denote the corresponding 
normalized, real-valued, and mutually orthogonal eigen- 
vectors. Let ipi = {ipii, ipi2, . • . , ipiw)'^ and ipi be the ith 
column vector of matrix 4'. Evidently, ^E* is an orthogonal 
matrix, satisfying 



(3) 



that is. 



N N f 1 if i = ■ 

y: ^..V',. = E = 0] ithJrwisa 

k—1 k—1 \ 

According to the properties of real symmetric matrices, 
one has 

*^F* = diag[Ai,A2,...,Ajv], (5) 

which, together with Eq. (|3]), leads to 

F = «'diag[Ai,A2,...,Ajv]«'^. (6) 

Equation ^ means that the entry Ty of matrix F has 
the following spectral form: 



N 

k=l 



Xklpkitpkj- 



(7) 



We now express matrix P in terms of and S. From 
Eqs. (dl and ^ one obtains 

P^S-^TS^ = S'-3<5diag[Ai,A2,...,AAr]*^S'5. (8) 

Let P* be the tth power of matrix P, whose ijth entry 
denoted by {p*')ij represents the probability for a walker 
to start from node i and reach node j in t steps. From 
Eq. ([U, one has 



F = S'-2i^S'-2 = S-^PS- 



(1) 



p'^s-H^\lMl\s^^ 



(9) 
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and 



By definition, the MFPT Fij can be evaluated as 



N 



fc=i 



In the limit, one has 



lim = i'liipijJ— = — 



(10) 



(11) 



where Eq. ([2]) and the relation w — ipi were used. There- 
fore, the stationary distribution for a random walk on 
weighted networks is 

TT = (7ri,7r2, . . . ,7rAr) =(—,—,..., — , (12) 

V s s s / 

which obviously depends on the strengths of the nodes 
and thus the weights of the edges. 



III. MFPT FROM ONE NODE TO ANOTHER 

MFPT is one of the most important quantities related 
to random walks, since it contains useful information 
about the random- walk dynamics 0. In what follows, 
we will derive MFPT for a random walk from one node 
to another in a weighted network, by using the formalism 
of generating functions ^47|] • 

Let {q*)ij represent the probability for a walker start- 
ing from node i to reach node j for the first time precisely 
in t steps. This first passage probability (q*)y relates to 
the transition probability (p*)ij by 



(13) 



Let Fij denote the MFPT for a random walk on a 
weighted network from node i to j. Then, 

oo 

F.j^Y.^iq'h. (14) 
t=l 

To find Fij, we introduce the following generating func- 
tions: 



QtA^) = XI (* 



t=0 



and 



P,,(x)=^(p*), 



(15) 



(16) 



t=o 



where \x\ < 1. The above two equations together yield 



Qij (*^) 



(17) 



d 



^3 - dx^^^^""^ 



(18) 



Thus, our goal is reduced to evaluating Qij{x) and then 
differentiating it. 

Plugging Eq. ^ into Eq. (HH) yields 



oo N 



Pij{x) = y^y^(Afcx)*V'fciV'fcjw — 

t=0 k=l ' 
AT — - AT oo 

= '^i'ki'4'kjJ — + y^y^(Afca;)*V'fci-0fcj./ — 



fe=l 



fc=l t=l 

oo N 



t=l fc=2 
N 



E(A 

.t=l 



^JT- ^ -; \ '^ki'ipkjy 

I — X ^-^ 1 — XkX 

k=2 



(19) 



where we have used Eqs. (j4|) and (|12p . In a similar way, 
Pjj (x) can be derived as follows: 



N 



t=0 fc=l 
N 



N 



^(Aix)x.-t-^^(Afcx)v; 



fc=l t=l 



N 



fe=i t=i 
1 ^ 1 



fc=2 Lt=l 



rk, 

(20) 



fc=2 



Substituting Eqs. ([H]) and (gHl) into Eq. dTT]) and carry- 
ing out some computation, we obtain 



N 



+ (1 - 2;) E 1 _ \ ^ "i^f^-^kjJ^ 
?..(2:) = ^^^^^ . (21) 



^ 1 

-. + (i--)Et3i^^: 



fc=2 



Differentiating Qij{x) with respect to a; and then setting 
x = 1, we obtain a closed- form expression for Fij as 



(22) 



which is useful for the following derivation. Similarly to 
the stationary distribution, the MFPT Fij is also related 
to the strength of the target node and weights of the 
edges involved. 
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IV. MFPT TO A GIVEN TARGET 



and 



N 



N 



After obtaining the MFPT Fij ft'om one node i to an- 
other one j in a generic weighted network G, we proceed 
to consider the trapping problem defined on G, which 
is a particular random walk with a perfect trap placed 
at a given node, e.g., node j. We use (Fj) to repre- 
sent the average trapping time, which is the mean of the 
MFPT Fij from a starting point i to the trap node j, 
with the starting point being taken over the stationary 
distribution [l^ that depends on the strength, as can 
be seen from Eq. ([T2l) . Note that the selection of the 
starting point is slightly different from that in most exist- 
ing works [lol - [l7j , where the distribution of the starting 
points is uniform. The choice of the stationary distribu- 
tion for the starting points adopted here is reasonable, 
since we focus on the trapping problem on weighted net- 
works. 

By definition, the ATT (Fj) is given by ^ 



1 



N 

1 i=l 



(23) 



Taking into account the expression of Fij in Eq. 
one straightforwardly arrives at 



N N 




Using Eq. 
Thus, Eq. 



O gives Ej^iV'fci 
24l) reduces to 



ki 



1 



^ 1 



fe=2 



kj 



(25) 



Equation ([25]) is a general result valid for the trapping 
problem on all weighted networks. 

We next derive a lower bound for (Fj). By Cauthy's 
inequality, one has 




N 



k=2 



\k=2 



(26) 

According to Eq. ([T]), every diagonal entry entry Sjj 



of matrix S satisfies 



Y.k=l ^k^lj 



> , J 



1,2,...,A^. Considering that the traces of matrices S 
and P are equivalent to each other, both equaling zero, 



we have s,- 



N 



N 



Ed 

fe=2 



0. Then, applying Ai = 1, one gets 

N 

\k)i>% = 1 - E ^^^fej- = 1 

(27) 



kj 



= E(i 

k=l 



k=l 



E ^fcj = E ^fej -^3^^- ^3 



(28) 



fc=2 fc=l 

Combining Eqs. (|25|) - ([28|) produces the following result: 



1 - TTj Sj 



(l-^,)^ = -(l-7r,) = --l. (29) 



3/ - ^ ~~\-^^"3J -—\-^~"J) 



This lower bound is strongly affected by the strength of 
the trapping node and the weights of all edges in the 
network. In addition, as will be shown in the following 
section, the lower bound for ATT can be achieved in some 
graphs. 



V. RESULTS FOR UNCORRELATED 
WEIGHTED NETWORKS 

We have shown that the primary quantities related to 
random walks on weighted networks rely on the node 
strengths and the weights of edges. Below, we will apply 
the above-obtained results to uncorrelated networks in 
order to uncover how these primary quantities change 
with the weight parameter 6. 

For uncorrelated networks, many quantities of interest 
can be determined explicitly. First, for the strength Si of 
node i, by definition we have 

= E id^d,f = {d,f E d,P{d'\d,){dy , (30) 

where P{d'\di) is the conditional probability |48j| that a 
node of degree di has a neighboring node with degree 
d'', dmin and dmax are the minimum and maximum node 
degrees, respectively. Let P{d) be the degree distribu- 
tion of the network. Then, for an uncorrelated network, 
Pid'\di) = d'P{d')/{d) and 

'f^^ {dr+^P{d') (d,)^+i(d^+i) 



id) 



id) 



(31) 

where {d^^^) is the {6 + l)th order moment of the degree 
distribution. 

The total strength s of all nodes can be evaluated as 

s = NP{s,)s, = NP{d,)s, = \, ^ , (32) 

(d) 

where P{si) is the distribution of node strengths that is 
equal to the degree distribution of the network. Thus, 
the stationary distribution for a random walk on an un- 
correlated weighted network becomes 



/Si 32 SAr\^ 
V S S 3 / 

^ f {d,Y+^ id2y+' 



{dN) 



e+1 • 



(33) 
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Equation (j33|) shows that the stationary distribution 
is dominated by the weight parameter 9. When 6* = 0, 
it reduces to the case of unbiased random walks. For the 
case oi 9 = —1, the stationary distribution is uniform 
with TT,; = 1/iV for node i. In the case of 6* > —1, tt^ 
increases with di, which means that it is easier to find 
the walker after a long time at a large-degree node than 
at a small-degree node. On the contrary, in the case 
oi 9 < —1, TTi is a decreasing function of di, and thus 
the final occupation probability for a node with a small 
degree is higher than another node with a large degree. 

From Eq. (1331) , we can also obtain the mean first return 
time (MFRT) Ru for uncorrelated weighted networks. 
By definition, MFRT Ru is the expected time for a walker 
originating from node i and returning to i for the first 
time. According to the Kac formula [M IH, [s^, the 
MFRT for a node coincides with the inverse probability 
to find the walker at this node in the final equilibrium 
state of the random- walk process. Thus, 



Ru = — 



N{d'^ 



e+i 



(34) 



which also depends on 9 and has been derived in 51] by 
using a different approach. 

Finally, for the lower bound of ATT on an uncorrelated 
network, inserting Eqs. 1^ and (HH) into Eq. yields 



> 



N{d^ 



{d,Y 



1. 



(35) 



This lower bound is sharp. It can be reached in an iV- 
node complete graph, where the MFPT from any node 
to j equals — 1 [53] , which is consistent with Eq. (1551) . 
Moreover, the ATT for complete graphs is independent 
of parameter 9. 

Different from complete graphs and some other regular 
graphs belonging to circular graphs, the lower bound for 
ATT provided by Eq. on non-regular networks relies 
heavily on the parameter 9. For those networks with 
small relaxation time when the stationary distribution is 
reached, the lower bound can be attained (see [53| for 
explanation) but the results for different destinations are 
closely related to 9. For 9 > the ATT for large- 
degree nodes is small, since Fj decreases with growing 
dj. For 9 = —1, the ATT is identical for all destination 
nodes irrespective of their degrees, while for 9 < —1, 
the walker can find small-degree nodes more easily than 
finding large-degree nodes, since in this case the walker 
is biased toward those nodes with small degrees. 

Note that for the 9 — Q case, it is reduced to unbi- 
ased random walks. For this particular case, Eq. (p9)) 
shows that the lower bound of ATT is proportional to 
the inverse degree of the trap node. Such a scaling has 
been reported in [31 and 53|. In [3], a main concern is 
ATT on scale-free networks, which are ubiquitous in real 
systems [s^ Issj. 

Next, we show that our results for ATT given in 
Eq. can display various different scalings when the 



trap is located on a target node having a maximum de- 
gree on a scale-free network and provide important in- 
formation about random walks on the scale-free network, 
encompassing that in [is'l as a particular case. 

For an uncorrelated scale-free network with a constant 
average degree {d) and a power-law degree distribution 
P{d) - d"^, the term (d^+i) in Eq. ^ can be evaluated 
as follows: for 9 = -1, (d^+i) = 1; for = 0, (d^+i) = {d) 
is exactly the average node degree; while for other cases 
(6* , -1), by integrating d^+'^P{d) from dmin to rfmax 
with respect to d, we obtain {d^'^^) to be approximatively 



(^max) 
((iniin) 



'-7+2 
61-7+2 



= 7-2, 
>7-2, 
< 7 - 2. 



(36) 



On the other hand, it has been established 42] that 
for an uncorrelated scale-free network with 7 > 3, its 
diaax approximately obeys the relation dmax A^i/(t~i). 
Thus, when a trap is fixed at a hub node having degree 
dmax on such a network with 7 > 3, the lower bound of 
ATT, (Fd^^^), can be expressed in terms of the network 
size N. For 9 = —1, (Fd^^^) ^ N — 1, which scales 
linearly with N. In fact, for this special case, the ATT is 
independent of the network type and the trap position. 
For 9 = 0, {Fd^^J = iV(T-2)/(T-i)^ growing sublinearly 
with N as observed in [l8!|. For 9 ^ , —1, we distinguish 
three cases: 



In N, 



= 7-2, 
>7-2, 
< 7 - 2. 



(37) 



Note that the approximate formulas given in Eq. (j37p are 
valid for small 9. Particularly, when 9 = 0, Eq. (1571) is in 
agreement with the previous result obtained in 18|. For 
other 9, Eq. (I37p may yield a large deviation from the 
true values [36]. 

To confirm our analytical results given in Eq. p7p . we 
perform extensive simulations of biased random walks 
on the Barabasi- Albert (B A) scale- free network |5^ , the 
degree distribution exponent 7 of which is equal to 3. 
In Fig. [1] we present our numerical results for the ATT 
to a most-connected node of BA networks with average 
node degree 4, which confirm the theoretical prediction 
provided by Eq. ([gT)). 

Equation (|37l) shows that, for large 9, one can find 
the trap hub easily. For example, in the whole range of 
— 1 < 6* < 0, when 9 is not very large, the ATT increases 
sublinearly with N, displaying a high trapping efficiency. 
Actually, for some particular 9 this result is also valid 
for correlated scale- free networks. Previous studies have 
indicated that, in the case oi 9 = 0, the ATT to a hub 
node in some correlated scale-free networks, such as the 
pseudofractal web [13, ill] , the Apollonian network [l^ , 
hierarchical scale- free networks [13, IH, [1^ , and modular 
scale- free networks [13, [13], behaves sublinearly with N. 
For small 9, especially for 9 < —1, Eq. (I37p shows that 
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FIG. 1: (Color online) Average trapping time to a hub node 
for biased random walks on the Barabasi- Albert scale-free net- 
work for various 9 and network size A*'. The symbols represent 
the numerical results generated by Eq. (|24p . while the solid 
lines stand for the corresponding analytical results given by 
Eq. dSZl). 



of random walks, the transition probability from node i 
to node j is proportional to the weight of the edge con- 
necting i and j, which is (didj)^ , where di and dj are, 
respectively, the degrees of i and j , and 6* is a control pa- 
rameter of weights. We have calculated analytically the 
expressions of the stationary distribution, MFPT from an 
arbitrary node to any other node, and ATT to a given 
trap node with a lower bound obtained. The resultant 
formulas for all these interesting quantities are sensitive 
to the change of 9, implying that 6 plays an essential 
role for random walks on weighted networks. For the 
special case of uncorrelated networks, we have derived a 
succinct expression for the stationary distribution, based 
on which we have further shown how the lower bound 
of ATT scales with the network size. Particularly, we 
have provided a comprehensive analysis for ATT on un- 
correlated scale-free networks, when a trap is positioned 
at a hub node. Our work provides an efficient method 
for controlling random- walk dynamics; using this method 
one can accelerate or decelerate the diffusion process on 
a weighted network. 
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VI. CONCLUSIONS 

In this paper, we have presented a unified framework 
for random walks on weighted networks. In the process 
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